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Abstract 

A relativistic model of a heat conducting collapsing star, which in- 
cludes thermal pre-relaxation processes, is presented. Particular attention 
is paid to the influence of a given parameter defined in terms of thermo- 
dynamic variables, on the outcome of evolution. Evaluation of the system 
when passing through a critical value of the aforesaid parameter, does not 
yield evidence of anomalous behaviour. 
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1 Introduction 



The behaviour of dissipative systems at the very moment when they 
depart from hydrostatic equilibrium has been recently studied [jj], ^, 

In |1| it appears that a parameter a formed by a specific combina- 
tion of thermal conductivity coefficient k, relaxation time t, temper- 
ature T, proper energy density p and pressure p, 

kT 



may critically affect the evolution of the object. Specifically, it was 
shown that in the equation of motion of any fluid element, the inertial 
mass density term vanishes for a = 1 (critical point) and is negative 
beyond that value. 

In some cases (pure shear or bulk viscosity) [Q, the critical point is 
well beyond the point where the causality requirements are violated 
and therefore forbidden. 

In others (pure thermal conduction) [jj], Q , later requirements are 
violated slightly below the critical point. 

In the general case ( heat conduction plus viscosity) it appears 
that causality may break down beyond the critical point [Q. 

However, it should be stressed that the critical point (as well as 
causality conditions) is obtained in the context of a linear perturba- 
tive scheme, where the system is evaluated immediatly after leaving 
the equilibrium and therefore time derivative of radial velocity as well 
as dissipative variables are considered small quantities, such that only 
linear expressions of them are kept. On the other hand the vanish- 
ing, at the critical point, of the inertial mass density term, leading to 
accelerated fluid elements in the absence of total radial forces, might 
suggest that linear approximation is not longer reliable there [^, g|. 

The question arises then, whether a physical system may reach the 
critical point without exhibiting a clear unphysical behavior or if, on 
the contrary, for any physical system reaching that point the march 
of physical variables becomes physically unacceptable as suggested 
by the linear approximation. In other words we want to see how the 
physical meaning of the above mentioned parameter a, as implied by 
the linear approximation, carries over to nonlinear regimes. 

In order to elucidate this question we shall consider here an exact 
numerical model of an evolving dissipative star (without viscosity). 
The system is forced to evolve through the critical point, and funda- 
mental variables are monitored to detect any anomalous behaviour. 
Modeling is performed by means of the HJR method [3], a brief re- 
sume of which is given in next section . 
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The heat conduction equation is given in section 3, and the model 
is described in section 4. Finally, a brief analysis of our results is 
given in the last section. 



2 The HJR method 

We shall consider a non-static spherically symmetric distribution of 
matter which consists of fluid, which may be locally anisotropic, and 
heat flow (radiation). Assuming Bondi coordinates [||, |6|, the metric 
takes the form 

ds 2 = e 2fs 



V 2 

— du + 2dudr 
r 



{.m- ...... 2 



+ sin 2 9d<p 2 ) , (1) 



where u is a time like coordinate , r is the null coordinate and 9 
and <f> are the usual angle coordinates. A generalization of the "mass 
aspect" defined by Bondi et al ||] can be introduced by means of 
function m(u, r) 

V = e 2/3 (r -2m(u,r)), (2) 

where /3 and V are functions of u and r. Inside the fluid distribution, 
the stress-energy tensor can be written as -see for details- 

Tfw = (p + Px)U„U v - Pxg„ v + (Pr - P±)x»X» + 1%nU v) , (3) 

where p, P r , P± are the energy density, radial pressure and tangen- 
tial pressure respectively as measured by a Minkowskian observer 
in the Lagrangean frame, and \n = —Qui It being q = ^J—q^q^ the 
heat flow. Using a Lorentz transformation in the radial direction and 
the coordinate transformation between Bondi coordinates and local 
Minkowskian coordinates 

dt = eP(^du+^dr), (4) 

dx = e P ^dr, (5) 

dy = rd9, (6) 

dz = rsinOdcj), (7) 

it is possible to express the components of the stress-energy tensor in 
Bondi coordinates, in terms of variables measured in the Minkowskian 
Lagrangean frame. In Bondi coordinates the four- velocity U^, and the 
heat flux vector q^ are given by 

1/2 

^ = e P U/T „ ;. a M/ 2 ,^(fr^) .0,0), (8) 
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and 

«"=*- i-^i^r^j^-^- < 9 > 

where co is the velocity of matter as measured by the Minkowskian 
observer defined by (Q)-(|7]). 

At the outside of the fluid distribution the space-time is described 
by the Vaidya metric which in Bondi coordinates is given by (3 — 
and V = r — 2m(u). 

The Einstein field equations, inside the matter distribution, can 
be written as 



(-m e~ 2/J + (1 - 2rh/r)m 1 ) = ^{p + 2uq + P r u 2 ), (10) 



47rr(r — 2m) 1 — uj 2 

fhi 
Airr 2 



I") 



doie~ its I, rh^, 9 #u 3/3i(l - 2mi) - TOn , N 

-—a + o-(l-2- (2/?ii+4/3?- — +— — ^ -=Pl, 13 

where subscripts and 1 denote partial derivative with respect to u 
and r respectively, and the effective energy density 

~p=-J—(p- q (l-u)-P rW ), (14) 

and the effective pressure 

P = T ^(-wp-g(l-w) + P r ) (15) 

are two auxiliary functions introduced in the HJR formalism |l(J 
whose physical meaning becomes clear in the static case, in which they 
reduce to the energy density and radial pressure respectively. 

Matching the Vaidya metric to the Bondi metric at the surface 
(r = a) of the fluid distribution implies see 
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for details 



Pa = -UaPai (16) 

which is equivalent to the well-known condition 

q a = Pra, (17) 



for radiative spheres [12|, (subscript a indicates that the quantity is 
evaluated at the boundary surface). 
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The HJR method is based on a system of three ordinary differential 
equations for quantities evaluated on the boundary surface (surface 
equations), which will allow us to find the evolution of the physical 
quantities. 

The dimensionless variables 



.4 



m(0) 



M 



m 
m(0) 



m(0) 



F= 1 



2M 



1 



1 



(18) 



are defined to derive the surface equations, where m(0) is the initial 
total mass of the system. Then, using ( [j"8| ) and boundary conditions 
one obtains the first surface equation 



A = F{n - i), 



(19) 



which gives the evolution of the radius of the star, (where dot denotes 
derivative with respect to the dimensionless u). 
The second surface equation 



F 



2L + F(l-F)(n-1) 
A : 



(20) 



gives the evolution of the redshift at the surface -see ^j, [n]] for 
details. The luminosity, L, as measured by an observer at rest at 
infinity reads 



L = -M = 



E 



= EF = E(2fl — 1)F = 4irA 2 q a (2f2 — 1) F, (21) 



where z a refers to the boundary gravitational redshift, E is the lu- 
minosity as seen by a comoving observer, and E is the luminosity 
measured by a non comoving observer located on the surface. 

The third surface equation is model dependent. For anisotropic 



fluids the relationship (T/f 



can be written as 



F U p a Pa 



-rFCl- 

A 



(i-«) 



AirApa 



being 



Pi. 



m - i 
n 



p + p 

1 - 2m/r 



F 



2A 



Pa 



Pa 
Ap a 



-Pr)a 



4 nr p + _ -(P x -P, 



(22) 



(23) 



Expression (22) generalizes the Tolman-Oppenheimer-Volkov equa- 
tion to the non-static radiative anisotropic case. 
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The HJR method [Q allows us to find non static solutions of the 
Einstein equations from static ones. The algorithm, extended for 



anisotropic fluids, can be found in [10 . Nevertheless, for our purpose 
here, it is only necessary to find the evolution of A, F, f2, and L. 
Thus, the algorithm can be resumed as follows 

1. Take a static but otherwise arbitrary interior solution of the 
Einstein equations for a spherically symmetric fluid distribution 
("seed" solution) 



P st = P(r), Pst = p(r). 



(24) 



Assume that the r-dependence of the effective quantities is the 
same as that of the energy density and radial pressure of the 
"seed" solution. Nevertheless, note that junction conditions in 
terms of effective variables, read as (16). This condition allows 
us to find out the relation between the u-dependence of p = p(u, r) 
and P = P(u, r). 



3. If we have an expression for p( 



u.r = a) 



and P(u, 



possible to solve the system of surface equations (19), (20), and 
(22) for a given luminosity. 



a), it is 



3 Heat conduction equation 

As mentioned in the Introduction, the value of parameter a, defined 
as 

(25) 



r(p + p) 

appears to be important in the evolution of radiating stars. In fact, 
as it is shown in [[]], W, if the system reaches the critical point (a = 1), 
then the inertial mass term vanishes. Furthermore it can be shown 



that conditions ensuring stability and causality [13, [L4j] are violated at 
the critical point [2]. This strange result may have, at least, two differ- 
ent interpretations. On one hand, it may happen that, because of the 
fact that the inertial mass term vanishes at the critical point, a linear 
perturbative scheme fails at the critical point, invalidating thereby 
any restriction obtained from linear aproximation (e.g. causality and 
stability conditions mentioned above). Alternatively, it may happen 
that the system is actually prevented from reaching the critical point, 
in which case it should exhibit some abnormal behaviour when ap- 
proaching the critical point. The correct interpretation can be found 
by solving the Einstein field equations, together with the transport 
equation, without applying linear perturbation theory. 



G 



In order to study the behaviour of the system at the critical point, 
it is convenient to adopt a physical framework in which condition 
a = 1 could be overtaken for reasonable values of physical quantities. 
A good candidate seems to be a collapsing neutron star in which 
the equilibrium is reached by means of a huge emission of neutrinos. 
The large temperature reached during the process of collapse would 
explain values of a greater than 1. 

We assume the evolution of the heat flow to be governed by the 
Maxwell-Cattaneo transport equation, 



= Kh} 1 



V {T,u 



-TU V 



(26) 



where 



U„= U a U v , a , (27) 
%= U a q,, a , (28) 

and k and r denote the thermal conductivity coefficient and relax- 
ation time respectively. Evaluating (^6|) in the surface, the transport 
equation can be written as 



rq a +q a ^F(2n-l) 



T a ~T la F(2n - 1) - T a 



Q(l - F) 
2A 



L 

AF 



n 



2il-l 



(29) 



The thermal conductivity coefficient for a mixture of matter and 
radiation is given by | 15 | 



K = -bT 3 r col , 



(30) 



where b = 7N v a/8 for neutrinos, being a the radiation constant, N v 
the number of neutrino flavors and r co / the matter-neutrinos time 
collision. 

On the other hand, the luminosity perceived by a comoving ob- 
server located on the surface can be connected with the effective 
temperature, T e ff, as 



L 



E = — = AnA 2 q a = Uir^aTtf] 



(31) 



where a = a/4 is the Steffan-Boltzman radiation constant. The con- 
cept of effective temperature has been widely used in theory of stellar 
atmospheres Jl^, p. 586], [0, p. 70], and Jl8|, p. 295]. The effective tem- 
perature can be connected to the material temperature by means of 



T 4 = - FT 4 1 

a r, L J -e//J r=a 



(32) 
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Thus, from ( pl| ) and (32), the dimensionless heat flow at the surface 
is 



la 



rpl 

i a' 



(33) 



29.- 1" 

where C = crM§ ~ 3.4097 x 1CT 54 K 4 , and £ = M /M Q . In the HJR 
formalism the initial mass is normalized to unity. Thus, every term 
in (E9H is dimensionless and k must be expressed as 
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(34) 



On the other hand, evaluating (E5J) at the surface, and by means 
of the definition of dl4|), ( fl7| ) and (|34|), it takes the form 

14y 4 $ 2 iV„n(2n - 1) , , 

Q 3(p a (2fi-l) 2 +4y 4 e 2 r!)' 1 j 

where we have assumed r ~ r co /, and the parameter y has been denned 
as 



y A = CT a 4 . 

Now, it is possible to write y in terms of a by means of ( |35| ) 

3p (2Q-l) 2 a 



2/ 



2^ 2 n(7JV w (2n-l)-6a) : 



and by substitution of this expression and (36), into ( |33| ) 

3p Q (2fi - l)a 



la 



0(7.^(20-1) -6a)' 



(36) 



(37) 



(38) 



After some elementary algebra, expression (29) can be written in 
terms of a instead of y and q a 



Hi 



F(2n _ 1) = _0_(_L_ 

y v ; 20- 1 V$ + 6a 



n(F-i) x 



3-y/F 



2A 



7t co1 N^2VL - 1 




1 



6a 



($ - 12aO) 
0$(20- 1) 



$ + 6a 



(39) 



where $ = 7N V (29 - 1) - 6a. 

Our purpose is to discern what anomalous effects (if any) can take 
place when the system overtakes the critical point. Thus, it seems 
reasonable to impose a profile for a, and study the evolution of the 
system. The Maxwell-Cattaneo transport equation (39) allows us 
to find the temperature gradient in the surface for a given a(u). At 
first glance, expression (39) does not seems to present anomalous 
behaviour in the critical point. Nevertheless, it is necessary to study 
the complete evolution of the system to confirm this suspicion. The 
method that we shall use can be summarized as follows 
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1. Impose a profile for a. According to (^5|), a must vanish initially 
if the system departs from equilibrium (T a (u = 0) ~ 0) 

2. For a given a, it is possible to find the luminosity L by means of 



(31) and (38) 



Follow the HJR method, outlined above, to solve the system of 
surface equations. 

Once A, F, and £1 are found it is possible to find the temperature 
gradient by means of (|39|). Assuming that cooling by absortion 
and emission of neutrinos is the responsible to drive the star 
to a new equilibrium state, the mean collision time r co / can be 
roughly expressed as [Q 

p\/yy 

where A = 10 9 K 3 / 2 m _1 , and y stands for the electron fraction. 



4 The model 

In order to study the evolution of the system beyond the critical point 



we adopt the Gokhroo &; Mehra-type solution [19| as the "seed" solu- 
tion. This solution corresponds to an anisotropic fluid with inhomo- 
geneous energy density, and it can take account of the large density 
and pressure gradients close to surface [Q. 

It can be shown [Q that the effective energy density and effective 
pressure are given for this model by 

P = *^(l (41) 



K V v 'A\ 
and 

where n > 1, and the central energy density in the static case p c is 
given by 

15 

Pc - 4*4,(5-3*,) (43) 

If we define 

1 = (44) 



then functions K(u), and G{u) are 



K{u) = < 



1 + Wl- 



( HK a \ (1 
{ 5 7 ) K 



1-W1 



and 



G(w) = 1 - 



- (^) (W 

(i-n)(i-J0 



if K > I 



ifK < 



(45) 



Fnx 



(46) 



Thus, the system of surface equations for this model is 

a= F(n-i), 



F=±[2L + F(l-F)(n-1)], 



4/clo 2 



F K (l-K) 3jKA 3 (2Cl- 1)(1-K) 



where 



A = - X) 



2K n 



30- 1\ 3 + F 2F0 



o 



2A A(l-K 



(47) 
(48) 

0(1 - 0)A, (49) 
(50) 



* = — ^A 7 A 2 X 2 (1 - G) n 
10K o 1 y ' 



A 2 
' 2F 



3 ^A 2 F 2 (1-G) 2 "-^F 2 (1-GV 



-i , 7 



4 1- 



32^ X(l-ir) 
K 

(51) 

The mean collision time (|40|) can be rewritten in terms of 0, a and 

p as 



T C ol = 



o$ 



y{$ + 3a) V Pa 



11/8 



20-1 



7/4 



2( 
3a 



-i 3/8 



whereas the luminosity can be found from (31) and (38) 



L : 



127rA 2 p a F(20 - 1) 
0$ 



(52) 



(53) 



Thus, the temperature gradient in the surface can be found by means 
of (52), ( |53| ) and ( |39| ) for a given N v . 

The system of equations (47-49) and (39) has been solved for two 
set of initial values: 
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1. Model with initial values A a = A(u = 0) = 6, il a = Sl(u = 0) = 1, 
n = 1, K = 0.999, A = 1/3, N v = 3 and y = 0.2. This initial 
configuration corresponds to a star initially at rest, with a radius 
of 8862 meters, a central density equal to 1.70 x 10 15 g cm -3 , and 
a surface density equal to 1.70 x 10 12 g cm -3 . Its initial mass is 
1M©, and its surface redshift is close to 0.225. We assume the 
center of the star composed by a highly relativistic Fermi gas. 
Thus, A = 1/3. The profile for a has been assumed as a gaussian 
centered in u = u pea k, i.e. 



exp 



tpeak 



(54) 



We have taken u peak = 100 (~ 0.49 msec) and A = 13 (~ 0.064 
msec). For this values the temperature changes considerably 
along 1 msec and a(0) ~ as is demanded by condition T(0) ~ 0. 
We have considered three profiles with a max = 0.9, 1 and 1.1 
(figure 1). 

2. Model initially in slow contraction with A a = 5 (~ 7838 meters), 
n o = 0.99926 (M < 1), ^ = n = 1, K a = 0.999, A = 1/3, N v = 3 and 
y = 0.3. These values correspond to p c ~ 2.94 x 10 15 g cm -3 and 
p a ~ 2.94 x 10 12 g cm -3 . The profile for a has been assumed as 



a = 2 — cxp 

and 



(55) 



a = 1. (56) 
As in the previous case A = 13 (<~ 0.064 msec) (figure 2). 



5 Discussion 

As mentioned before, our purpose here is to clarify the physical sig- 
nificance (if any) of the critical point in the evolution of a dissipative 
self-gravitating system. 

In particular we want to find out if at the critical point the system, 
described by the full theory (without any approximation), exhibits 
an anomalous behaviour, as might suggest a linear approximation 
approach. 

The three profiles of a considered are displayed in figure 1. In one 
case the system never reaches the critical point, in other case it is at 
the critical point at some moment of its evolution, and in the third 
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case the system goes beyond the critical point before returning to 
equilibrium. 

Figures (3)- (6) show the evolution of temperature, temperature 
gradient, radius and surface velocity for the three different profiles of 
a. 

It is apparent from these figures that even though increasing values 
of a are associated with more unstable configurations (in the sense 
of faster collapse), nothing strange seems to happens at or beyond 
the critical point. Figure (7), showing the evolution of the ratio of 
the neutrino mean free path to the radius of the star, indicates that 
the diffusion approximation is valid during most part of the emission 
process. 

To reinforce this conclusion we have performed another numerical 
simulation with the profiles of a given in figure (2). The evolution 
of relevant physical variables displayed in figures (8)-(ll), confirms 
the conclusion emerging from figures (3)- (6). The validity of diffusion 
approximation is corroborated in fig. (12). 

To conclude, we may say that the increasing of instability (in the 
sense mentioned above) with higher values of a, seems to confirm the 
decreasing of the inertial mass density factor, predicted in the linear 
approximation. 

However, nothing dramatic appears at the critical point in the 
exact modeling, indicating that the later approximation is not reliable 
at (or close to) the critical point, as suggested before |^]. 
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Figure captions 



Figure 1.- Profiles of a given by expression (54). The values of a max 
are 0.9, 1 and 1.1. 

Figure 2.- Profiles of a for the model initially in slow contraction ( |55| ) 
and (|56|). 

Figure 3.- Temperature evolution corresponding to three profiles of 
a shown in figure 1. 

Figure 4.- Temperature gradient for the model initially at rest. 

Figure 5.- Evolution of the radius for the model initially at rest. 

Figure 6.- Evolution of the surface velocity for the model initially at 
rest. 

Figure 7.- Mean free path of neutrinos in the first model. 

Figure 8.- Temperature evolution for the model initially in slow con- 
traction. The profiles of a corresponding to this model are shown 
in figure 2. The curve labeled with a = 1 corresponds to the pro- 
file (56), whereas the other one corresponds to the profile of a 
given by expression (|55|). 

Figure 9.- Temperature gradient for the model initially in slow con- 
traction. 

Figure 10.- Evolution of the radius for the model initially in slow 
contraction. 

Figure 11.- Evolution of the surface velocity for the model initially 
in slow contraction. 

Figure 12.- Same as figure 7 for the model initially in slow contrac- 
tion. 
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